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Abstract 

We propose a model to analyze the insurgence of pull-in and wrinkling failures in electroactive thin films. We take 
in consideration both cases of voltage and charge control, and study the role of prestretch and size of activated 
regions, which are essential in the analysis of realistic applications of EAPs. Based on simple geometrical and 
material assumptions we deduce an explicit analytical description of these phenomena, allowing a clear physical 
interpretation. Despite our simplifying assumptions, the comparison with experiments shows a satisfying qualita- 
tive and, interestingly, quantitative agreement. In particular our model shows, in accordance with experiments, the 
existence of different optimal prestretch values, depending on the choice of the actuating parameter of the EAP. 

1 Introduction 

Electro-active polymers (EAPs) are very promising materials for several technological applications in the 
robotic, medical, and biological fields [Hill]. These materials are characterized by important qualities 
such as lightweight, small size, low-cost, flexibility, fast response. The rapid technological advances in 
material science have recently improved EAPs performances in terms of their traditional limits, such as 
small actuation forces, low robustness and the requirement of very high electric fields. A prototypical 
device for electro-actuation based on EAPs is constituted by a thin layer of polymeric material sandwiched 
between two compliant electrodes. Application of a voltage to the two parallel electrodes generates an 
electric force (measured through the so called Maxwell stress [5]) which induces a compression of the 
layer; since the polymer is typically nearly incompressible, the resulting transversal extension is used as 
a mean of actuation. 

A typical drawback in technological applications based on the described electro-activated films is the 
insurgence of different types of instability, leading to failure (e.g. [15]). One of the well known phenomena 
occurring in these devices is the so called pull-in instability. In order to describe this type of instability, 
one may consider a toy model composed by two rigid conducting plates connected by an insulated linear 
spring, with a rest distance equal to d. If the plates are subjected to a voltage V, the Coulomb forces 
between the plates will tend to attract them so that the spring will result compressed. It is easy to check 
that there exists a voltage threshold V* such that for larger voltages no more equilibrium is possible 
between the elastic and electric forces and the top plate slams onto the bottom plate. The resulting 
thumb rule used in EAPs and similar MEMS devices is that the critical distance associated with V* 
equals d/3. 

Clearly, the quantitative determination of the pull-in voltage is not a simple task for real dcformablc 
devices: the failure mechanism of electro-activated thin films depends not only on electro-elastic interac- 
tions, but also on the actuator shape and on purely electrical effects. A further purely mechanical failure 
mechanism in EAPs is represented by the insurgence of compressive stresses inside the polymer film. In 
typical actuation systems (e.g. [1 1| . |15j ) the polymeric film is constrained inside a rigid frame which serves 
to assembly the actuator, to provide electrical insulation, to transmit the actuation force, and to provide 
a pre-stretching of the polymer. Under the action of electric forces the constrained polymer film may 
undergo in-plane compressions, typically shed in evidence by the insurgence of a buckling type instability 
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known as wrinkling in membrane theory ( [161 Another crucial aspect for technological applications 
is given by the observation that the polymer film is tipically activated only in an internal region, due 
again to both assembly and insulation motivations. 

In the recent past several theoretical and numerical analysis have been proposed in order to elucidate 
the causes of the various failure mechanisms of EAPs actuation devices. The importance of prestraining 
has been clearly evidenced in and |14j , where experimental observations have shown the existence of 
optimal prestraining values. An interesting analysis of the benefits of prestraining has been carried on 
in |10) for an unconstrained layer activated on the whole surface; in the same work, the necessity of a 
correct description of the strongly non linear material behavior of the polymer has also been evidenced. 
A similar problem is considered in [13] where the authors have discussed the relevance of an energetic 
analysis in the description of the pull-in instability. A numerical analysis for a circular planar actuator 
has been recently presented in |15| . where the main role of prestretch and stretch rate up to failure 
are numerically analyzed. An energetic analysis of prestress and pull-in for a homogeneously deformed 
dielectric elastomers is studied in [3D]. Finally, a theoretical analysis of the insurgence of deformation 
localization in a variational framework was recently proposed in |17| and |21] , where the role of damage 
and dissipation were considered. 

In this paper we propose a simple prototypical model for a slender prcstraincd actuator with an 
inner region activated by assigned voltage or charge (see the scheme in FigfT]). Based on simplifying 
geometrical and material hypotheses we elucidate the combined roles of prestrain, material properties 
and dimensions with respect to the wrinkling and pull-in failure modes. Due to the film thinness we 
neglect bending resistance and the formation of a boundary layer between activated and non activated 
regions, so that we assume a sharp interface between such regions [5] 0] and we ignore fringing effects 
[12| . Our geometrical assumptions are well suited for linear actuators typically adopted in artificial muscle 
technologies [14] . Of course our approach can be extended to more complex geometrical schemes and 
constitutive assumptions by using standards numerical methods. 
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Figure 1: Scheme of the model and boundary conditions. The gray region represents the activated region. 



The paper is organized as follows. In section [2] we recall the basic electromechanical equilibrium 
equations. In section [3] we introduce our prototypical example for a thin electroactivated system. In 
section [4] we obtain analytical solutions for both boundary value problems with assigned charge and 
voltage, as a function of given prestretch. Then, in section [5] we discuss the results of our model. In 
particular, we deduce analytically the existence of an optimal pre-strain (experimentally detected in 
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|14j ) which depends on the material modulus of the polymeric layer and on the size of the activated 
region. Shortly, the model describes how prestretch has a positive effect on the wrinkling instability 
and a negative effect on the pull-in instability, leading to the existence of an optimal prestretch. The 
corresponding value sof the activation parameters are then deduced. In section [5] we show that the 
model well describes the experimental behavior. Finally, in the Appendix we deduce a necessary stability 
condition for the electroactivc device here considered. 



2 Basic Equations 

In this section we collect the main equations that are at the base of our theoretical approach. We refer 
the reader to the recent paper [2] and to the references therein for details. Let X £ Bo and x € B be the 
typical point in the reference Bo and in the current configuration B of a continuous body, with x = f(X), 
where / is the deformation field. Moreover, we denote by (j>, D and E the potential, the displacement 
vector and the electric field in the current configuration, respectively. For a linear, homogeneous and 
isotropic dielectric we have 

D = eE = -egrad^, (2.1) 

where e = £o£d with £q the permittivity of free space and Ed the dielectric constant of the material. The 
equilibrium equations for the electromechanical problem are 

div.D = p, divT = b in B 

(2.2) 

[Z>] ■ n = q in dB q , Tn = s in dB s . 

Here p and q (assigned on dB q ) are the bulk and surface charge densities and b and s (assigned on 
dB s ) are the bulk and surface forces densities in the current configuration. T is the (symmetric) total 
stress tensor, that can be decomposed as the sum of the elastic (ela) term T ela and of the electric (elc) 
component T e ' c (i.e. the Maxwell stress) 

r|-i rj^ela _|_ rj-ie/c ^2 2) 

with 

T elc = e(E®E- \(E ■ E)I), 

(2.4) 

T$° = \E i E j -{E.E)6 ij , 

where i,j range from 1 to 3 (here and in the following we will use, when necessary, both component and 
absolute notations). The corresponding nominal (Piola-Kirchhoff) stress can be obtained as 

S = S e ' a + S e/C = jT e/Q F~ T + JT e ' c F~ T , (2.5) 
where F := V/ and J := detF. 



3 Setting of the problem 

We consider a dielectric elastomer actuator which occupies, in its natural configuration, a right prismatic 
region [0, £] x [0, L] x [0, H] . We use the lenght I of the elastomer to adimensionalize the thickness H = H/£ 
and the width L = L/l and thus we consider the prismatic region [0, 1] x [0, L] x [0, H] (see FigUJ where 
we also indicate the chosen coordinate system). In particular, we assume that the dielectric is thin, in 
the sense that H -C L <C 1. For thin films the membrane approximation of the nonlinear elasticity may 
be adopted. The standard membrane assumption (e.g. [3]) states that the bending stiffness is zero and 
any in-plane compressive stress leads immediately to the membrane buckling, with the appearance of a 
wrinkled regions. The presence of wrinkles is frequently encountered in experiments on thin dielectric 
elastomers and it represents one of the main observed failure mechanism of such devices. 
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We assume that voltage or charge can be controlled on a region of the thin film, more precisely on 
opposite strips of length \ 5; 1 belonging to the upper and lower faces of the actuator. This region is 
placed centrally in direction of e\ length and it covers the full length L in direction e2- The prismatic 
region of the actuator with upper and lower faces coinciding with the electrically controlled strips will 
be shortly denoted by active region, whereas the remaining part is simply called non active region. In 
accordance with the technological applications described in the introduction, we consider the possibility 
of a prestretch of the actuator along the e\ direction (see again FigfT]). 

We describe the layer deformation by two homogeneous deformations with eigenvectors coinciding 
with the coordinate axes. More specifically, we assume in the non active region a deformation / such 
that 

non active region: 



x\ = \\X\ 

x 2 = X2X2 , F := V/ = 
X3 = A3X3 



Ai 
A 2 
A 3 



(3.1) 



In the activated region we denote by X = (X\, A 2 , A3) and x = (xi,X2,X3) the generic point in the 
reference and current configuration, respectively, and by / the deformation. For this region we assume 



active region: 

x\ = \\X\ 
x 2 = A 2 A 2 , 
£3 = A3X3 



V/ 



Ai 
A 2 
A 3 



(3.2) 



The assumption of piecewise homogeneous deformation can be justified by the following hypothe- 
ses: piecewise homogeneous loading, thinness hypothesis, and constitutive assumption of homogeneous, 
isotropic, material behavior of the polymer. This type of deformations (which essentially neglects the 
presence of a geometrically-compatibile boundary layer at the interface between active and non active 
regions) is classical within the framework of phase-transitions [5] and of thin bodies subjected to homo- 
geneous boundary loading (see, e.g., [5j [4j ) . The balance equation at the interface may be formulated in 
the weak form [SJei = 0, e\ being the normal to the interface in the reference configuration. 

As well as the geometrical boundary layer we here neglect edge effects of the electric field, known in 
literature as fringing fields. Across the interface between the active and non active regions any disconti- 
nuity of the electric field \E\ must fulfill the jump condition (see e.g. [2]) e\ x \E\ = 0. The assumption 
of negligible fringing fields leads to an apparent violation of the jump condition, since it results that 
E = Ee$ in the active region whereas E = in the non active region, so that the component of E which 
is parallel to the interface is not continuous through the interface; this apparent violation is actually 
resolved by considering that in reality the electric field undergoes an abrupt, but continuous, variation 
from £ to in a region of negligible width. Also the assumption regarding electrical edge effects are 
acceptable, granted the thickness of the EAP is much smaller than its other dimensions (see e.g. [T2"]). 

Regarding the polymer material behavior we assume that this is incompressible, so that its deforma- 
tions respect the isochoric constraint 



A1A2A3 



1. 



AiA 2 A3 — 1. 



(3.3) 



Moreover, in a typical application scheme, the actuator is pre-stretched in the direction X\, say of an 
amount equal to v. By imposing that the total length of the actuator in the current configuration equals 
v, we get 

^ = xAi + (1-x)Al (3.4) 

To obtain analytic results we assume that the incompressible polymer is characterized by a neo-Hookcan 
constitutive response 

T ela = -ttI + fiB, S e ' a = T e ' a F~ T , 



ij 



ela 



-7fA- Sij + flXi, 



(3.5) 
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where 7r is a Lagrange multiplier due to the incomprcssibility constraint, fi is the shear modulus and 
B := FF T is the left 
stresses and pressure 



B := FF T is the left Cauchy-Green deformation tensor. By introducing the adimensionalized principal 



_ T\ 7[ 

n n n 

the non zero stress components in the non active region can be calculated by substitution of the defor- 
mation (|3.1[) in the constitutive response, which gives 

non active region: 

(3.6) 



s, 



= -7rAr 1 +A i , U = -tt + A? (i = 1,2,3). 



As anticipated above, we assume that the electric field in the active region coincides with E = (0, 0, E). 
By using (|2.3[) and (|2.4[) the principal stresses in the active region are given by 

active region: 



t« = -tt + A? +e£ (3-7) 



Si ' ~ ~ '% ft, 

The nature of the electric stresses, as well as their geometrical coupling with the deformation, changes 
according to the fact that voltage or electric charge are externally controlled by the device which drives 
the actuator. In the first case, the applied voltage V determines the electric field: 

eV V 
voltage control : D = (0, 0,— -), E =(0,0, —r-). (3.8) 

In the second case, the total electric charge Q is controlled on the electrodes, so that 

charge control :D = (0,0,^-2— ), E = (0,0,-J? — ). (3.9) 
AiA 2 ^4 eAiA 2 ^4 

Here A = x L is the reference area where the charge is applied. In the following the two driving mecha- 
nisms will be considered separately. 



4 Equilibrium for the voltage and charge controls 



We first consider the following pre-stretch deformation f p in the e\ direction 



vXi 
V2X2 
V3X3 



F p := V X / P = 



^00 
u 2 
is 3 



(4.1) 



where v = {v^v^} is the assigned longitudinal prestretch. (Here and in the following, for simplicity of 
notation, we sometime drop the index 1 when there is no ambiguity.) 

The principal Cauchy stresses in the prestretched configuration are given by 



iP - 



so that after imposing the conditions of vanishing stresses on the faces orthogonal to the axis X2 and X3 



t 1 ' 
'2 



0, 



(4.2) 



we obtain ir p = v 1 and v 2 
configuration amounts to 



"3 



v x / 2 . Thus, the only non-zero stress component in the prestretched 
p = v 2 -v-\ s\ = v-v- 2 . (4.3) 
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Now we consider the electric activation of the EAP. Regarding the nonactive region the boundary 
conditions t<i = £3 = give (Ai,A2,A3) = (A, A -1 / 2 , A -1 / 2 ) so that the stress in this region as in Eq.ns 
(|4.3[) is given by 

U = (A 2 - \- x )6a, Si = (A - A- 2 )^. (4.4) 
Concerning the electrically active region, after substitution of Eq.ns p.8[) in Eq. (|3.7p . we get 

1 eV 2 

voltage control: U = —tt + A 2 + — (26 i3 - 1), (4.5) 

2// X 2 H 2 

and, by using (|3.3|) . (13.71) and p.9[) we get 

1 Q 2 A 2 

charge control : U = -tt + A 2 + ^"^F^ia - !)• (4-6) 

In order that the EAP maintains its total length equal to the given prestretch also after the application 
of charge or voltage (see Eq. (|3.4[0 . we impose 

I/ = (1-X)A + X A. (4.7) 

This equation gives the following relation 

A=^, (4.8) 
1 - X 

so that the condition A > implies that 

A < v/ X (4.9) 



Using p~4|) and (|478|) we then obtain 

Let us now take in consideration the case when voltage is applied. In this case it is convenient to 
introduce the parameter k v defined by 

1 eV 2 



1 2ft H 2 ' 

which represents the (adimensionalized) electric free energy density under the assigned potential. By 
imposing the boundary conditions £2 = £3 = 0, after simple calculations we obtain the equilibrium stress 
and strains 

s = X- 



voltage control: < 



A 2 = ; = , (4.u) 

AVI - 2k v X 2 y ' 



\ 2 ^/1-2k v X 2 
1 



~ 2 _ sfl - 2k v X 2 

Ao — . 

A 

When the total charge is assigned, after imposing the boundary conditions £2 =£3 = and by introducing 
the (adimensionalized) electric energy density k q 

1 Q 2 
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wc obtain 



= A - 



s/1 + 2k q 



charge control: < 



2 _ VTT^ 

An — 



A 



A 
1 



(4.12) 



Ayl + 2k 

Finally, in order to deduce the equilibrium configurations we must impose the mechanical balance at 
the interface between the active and nonactivc regions, which reads 

s = s. (4.13) 

This equation let us determine A and, hence, the equilibrium configurations of the system. 



5 Wrinkling and pull-in instabilities 

The two failure mechanisms of wrinkling and pull-in will be considered separately for the cases of voltage 
and charge control. 



5.1 Voltage Control 

When voltage is assigned (see Figf2]) the equilibrium configurations, according with (|4.13[) . (|4.10l) . (|4.11[) 
correspond (if these exist) to the zeros of the function 

s " v ' 



qW = a — 



AV1-2k,,A 2 



measuring the stress discontinuity at the interface. In order to obtain physically reasonable deformations 
(see (|4.9p and (|4.11[l i ) we are interested to the solutions of g(X) = in the range (0, A ) where A is 
the min{l/-y/2K r , v/\}. We observe that the stress s in (|4.10[) is a (^-dependent) monotonic decreasing 
function of A, whereas an analysis of (|4.11j) shows that s is a (k v -dependent) concave function of A, 
defined for A G (0, l/y/2n v ), with s — > — oo at the boundaries of the definition domain. An inspection of 
the stress plots reported in FigJSJj shows that for given v and \ the condition g(X) = may have two 
roots A_ < A + (equal or distinct) or no roots in the range of interest. If two roots exist it is possible to 
show that only one of them is stable in a sense which is discussed in detail in the Appendix; by means of 
this condition, it easy to check that the stable solution corresponds to the root A_ . 

For given \, v we define Ky mp the value of the loading parameter in correspondence of which the stress 
evaluated at A_ is such that s = s = (see FigJH,): for values oi k v > n c ° mv the equilibrium stress is 
negative and, due to the lack of bending resistance, the film exhibits wrinkling, determining the device 
failure. Always keeping is, x fixed, a second eventuality corresponds to values of the load parameter for 
which the two solutions coalesce (see again FiglJ],), so that for values of k v larger than this critical value 
kI"" no more roots of s = s exist: this instance represents the pull-in instability. For n v > k p ^" the 
elastic response of the film is no more capable of balancing the compressive electric forces and the facing 
electrodes smash onto each other. Clearly, in correspondence of n v = k^'" it results g = and dg/d\ = 0. 

The situation is depicted in FigEl where both values of k c °"' p and of K"y" are plotted with respect to 
the prestretch v. Under the particularly simple assumptions of our model, the expression of k™"' p can be 
deduced analytically by imposing s\ = s\ = 0, which gives 

= \ (Ut - T±rX with %(v) = v -±*^±. 

2 \X 2 (v) X s (is) J X 
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(X = 0.5) 




Figure 2: Stress strain curves in the activated and non activated region in the case of assigned voltage: 
continuous lines represent the stress si(X,k v ) in (|4.11[) , of the activated region, whereas dashed lines 
represent the stress Si(A, v) in (|4.10[) in the non active region. In b) we show the attainment, for a fixed 
v of both the pull-in and wrinkling instability. In c) we show the stress curves si(A, v) corresponding to 
the maximum attained value of Ky 



The function k c °™ t (v) is plotted with a bold line in FigOJ it attains a maximum in correspondence of a 
prestretch 

vZ7 = 1 + (2 1/3 - 1)X- (5.1) 

The maximum value of n v which can be applied to the film without inducing compression then amounts 
to 

«5Z.= |(2~ 2 / 3 - 2- 8 / 3 ), (5.2) 

attained when 

A = A""" P = 2 1/3 . 

opt 

The analytical expression of K*f" is not easy to get, but the numerical solution of the conditions g = 0, 
dg/dX = gives k p J" for each value of v, which is represented by a light line in Fig|3l It is possible to 
show that, for fixed Xj the curve K"° mv (y) is always under the curve K p y"(u) but in a common point where 



8 



1 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 

Prestrctch v 



Figure 3: Critical values of the activation energy k v . Bold line represents the function k c ° mT {v) corre- 
sponding to wrinkling instability, whereas light line represent the function k\ ull {v) corresponding to pull-in 
instability. Dashed lines correspond to unstable configurations. 



these curves are tangent, corresponding to the simultaneous attainment of both wrinkling and pull-in 
failures. This point can be analytically determined by imposing that s = s = and dg/dX = 0, which 
after simple algebraic manipulations yields the value of the prestretch 

< = xiKKx) - 1) + 1, (5-3) 

where 

K1(X) ■■= (j^) 1 * ■ (5-4) 

It is easy to check that for all \ it results v"°™ v < v'" o f t l . 

It is interesting to observe that not all the values of the functions K^"(v) and K c ° mp (v) are actually 
attainable; regarding k^""(i/), the analysis of the functions s(A) and s(A) shows that for v > the 
pull-in failure is attained at positive stress, whereas for v < u^f it is attained at negative stress, which is 
not possible in our model (dashed light line in Figj3j). On the other side it is possible to check that the 
values of n c ° mp corresponding to v > v^" are relative to the unstable solutions A+ and hence these are 
not attainable (dashed bold line in Figj3|). As a result, the film failure is due to wrinkling for v < v^" 
and it is due to pull-in for v > v vu ^. As it is clear from the plot, the highest value of the load parameter 
which can be applied without generating failure of both types coincides with the maximum value of K c ° mv 
given in Eq l5.2[ which is attained in correspondence of v = ^™ f mp ; in this sense ^° t mp is optimal if one is 
interested at the maximization of the stored electric energy. 

A remarkable fact is that if one is interested at other useful activation parameters, such as the 
maximum actuating force or the maximum applicable electric field, then v"°™ v is no more optimal. This 
can be seen by evaluating the values of the stress gap with respect to the prestretched configuration, 
given by AS := n[s(v) — s(A)] and of the electric field E in correspondence of the initiation of both 
wrinkling and pull-in failures; as it is evident from the plots, the maximum values of both parameters 
are now attained for v = v r o u " tl that is for A = A^f (see Eqn.s (|5.4[) - ([5.3[) ). We remind that this value 
of v corresponds to the simultaneous attainment of the wrinkling and pull-in failures. After a short 
calculation, we deduce that the optimal value of AS 

AS-" = M<-(<)" 2 ), (5.5) 
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whereas the optimal value of the electric field is given by 

/ 3/i 



E C=\I— ((l-x) _ *(4 -*)-*). 
It should be underlined, in conclusion, that all the activation parameters AS, k v and E depend on 

(X = 0.5) 



AS p ""unstable 

A5P""stable 

AS c<,m, >unstable 




1.2 1.3 1.4 

Prestretch v 

(X = 0.5) 



E» u " unstable 
E pu " stable 
E ra "» stable 

B cmp unstab [ e 




1.4 1.6 1.; 

Prestretch 



a) 



b) 



Figure 4: Critical values of the activation stress AS a) and of the electric field E. Bold lines represent 
the values corresponding to wrinkling instability, whereas light lines represent the values corresponding 
to pull-in instability. Dashed lines correspond to unstable equilibrium states. 



the elastic modulus fj., and that AS and E grow with the fraction \ °f activated film. The resulting 
dependence of v pu " and i/ comp with respect to the fraction of activated material is summarized in Figj5] 
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Figure 5: Optimal prestretch (|5.1[) and (|5.3[) as functions of the size \ °f the activation region. 



5.2 Charge Control 

In the case of charge control, it is interesting to observe that the pull-in failure is not possible (this is 
in accordance with previous analyses [H]). This is due to the fact that for any v and \ the stress s 
is a decreasing function of A, whereas the stress s is an increasing function of A. On the contrary the 
wrinkling failure is possible and it is attained as soon as s(y, x) = s(k q ) = 0, that is 

a = K7» = i + — 

X 

which corresponds to a maximum load parameter 

^comp 1 ( /\comp\6 l\ 

K Q. m a x - 2 \\ A opt ) L j ■ 

Since the load parameter is an increasing function of v, it is clear that in the case of charge control it is 
not possible to find an optimal value of the prestretch. 

6 Comparison with experiments 

In this section we briefly compare the theoretical results of our model with the experimental behavior of 
electroactive polymer actuators 1 1 1 j . 

In FigfTk we reproduce the experimental behavior of silicone rubber films (Nusil CF19-2186) as re- 
ported in [11] , where circles represent the values of the electric field at failure for different assigned pre- 
stretches. In FigJT], we reproduce the results of our model, where we have considered an elastic modulus 
/i = lAIPa, as suggested by [IT] for this material. It should be remarked that the considered experiments 
refer to circular actuators, whereas we here have focused attention on linear actuators. Nevertheless, the 
comparison between experiments and theoretical results are in nice quantitative agreement. 

Quite interestingly, a careful analysis of the experimental FigJTb in light of the theoretical FigfTl, shows 
the existence of different behaviors, separated by a threshold prestretch which approximately corresponds 
to z/™ t ". For values of v lower than w 1.4 both the experimental and the theoretical results show an initial 
growth of the maximum applicable electric field with growing prestretch. According to our model this 
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(X = 0.5) 




0.5 1 1.5 2 2.5 



Figure 6: Stress strain curves in the activated and non activated region in the case of assigned charge: 
continuous lines represent the stress si(A,k ) in (|4.12|) . of the activated region, whereas dashed lines 
represent the stress si(A, v) in (I4.10P in the non active region. 



may be due to the initially positive effect of prestretch in avoiding the occurrence of wrinkling. Such 
positive effect, both in experiments and in our model, ceases to exist when the prestretch overcomes 
the critical value of about 1.4, since here the thickness reduction anticipates the occurrence of pull-in 
failure. A finer description of the film behavior for larger values of the prestretch certainly requires a 
better description of the exact film geometry and of the material response. To the price of renouncing to 
closed form solutions (which arc generally possible under very simple assumptions such those embraced in 
this article) , the methods we have discussed can be numerically applied to more complex geometries and 
material behaviors suitable for polymers undergoing large deformations, such as Mooney-Rivlin, Gent, 
Arruda-Boyce and others; some of these applications will be considered in a forthcoming paper. 



APPENDIX 



In this Appendix a necessary stability condition is obtained for the equilibrium configurations of a linear 
actuator undergoing piecewise homogeneous deformations. We observe that the homogeneity of the 
deformations in both the active and non active regions, together with the constraint of assigned prestretch 
v, essentially reduces the present equilibrium problem to a one-dimensional problem. Under conservative 
hypotheses, the total energy <§ stored in the device at given A can be obtained by integration of the 
internal working density along an arbitrary deformation path from the undeformcd configuration; since 
in this case the only non zero stress component is s we have 



X / s(x) dx + (1 — x) / s ( x ) dx 



HL 

This function must be minimized under the compatibility constraint v = xA + (1 — x)A, that gives 



A = A(A) 



i-x 
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Figure 7: Failure value of the electric field as a function of the prestretch: a) experiment (deduced from 
[TT]). b) theoretical. 



Let then A define a given configuration of the system, and let A + £ be a perturbed configuration that is 
obtained by shifting the interface between the active and non active regions. The perturbed energy is 

s'Cx + O f ' x+c - /" a ^ c t^ 

= X I s(x) dx + (1 — x) / s(x)dx. 



HL 

Thus, a stable configuration A must verify 

S(X + C) > <?(A) for all C 
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By considering the following expansion 



<T(A + C) 



<?(A) + Cx s(A)-s(A) 




2 _ ds(X) ds(X(X)) 
2 dX dX 



+ o(C 2 ), 



the vanishing of the first variation gives the Euler Lagrange equilibrium condition 



s(A) = s(A), 



whereas the positivity of the second variation gives 



rfi(A) ds{X{X)) 
dX ~ dX 



which is the necessary stability condition adopted in the paper. 
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